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Based  on  summarizing  the  latest  research  of  vehicle  exhaust  thermoelectric  generator  (TEG),  this  paper 
proposes  a  newly  designed  TEG  with  cylindrical  shell  and  straight  fins.  A  two-dimensional  heat  transfer 
numerical  model  using  finite  difference  method  is  established  under  steady-state  conditions.  This  model 
is  used  to  predict  the  output  performance  of  TEG  in  different  simulation  conditions.  The  average 
contribution  rate  to  the  output  power  CRvar  is  given  in  this  paper  to  evaluate  the  influence  level  of 
different  input  parameters.  Simulation  results  show  that  increasing  the  inlet  temperature  and  the  inlet 
mass  flow  rate  of  the  exhaust,  as  well  as  reducing  the  inlet  temperature  of  the  coolant  are  useful  methods 
to  improve  the  performance  of  TEG  in  a  certain  extent.  The  heat  transfer  numerical  model  presented  in 
this  paper  can  be  also  applied  to  the  numerical  simulations  and  the  optimization  work  on  early  designing 
stages  of  other  types  of  TEG. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

As  a  new  form  of  waste  heat  recovery  application,  thermoelec¬ 
tric  generator  (TEG)  can  absorb  heat  when  there  exists  temperature 
difference,  and  transfer  the  heat  directly  into  electricity,  with  the 
properties  of  no  pollution,  no  noise,  no  moving  parts,  high  reli¬ 
ability,  long  life,  etc.  This  technology  can  be  used  in  medical,  mili¬ 
tary,  energy  industry,  space  technology  and  other  fields  [1,2]. 
Vehicle  engine,  as  a  typical  representative  of  power  machinery,  has 
the  conversion  efficiency  of  fuel  combustion  only  about  40%,  and 
about  30%— 45%  of  the  total  energy  is  taken  away  by  engine  exhaust. 
When  a  TEG  is  installed  on  the  vehicle  exhaust  tube,  the  waste  heat 
can  be  recovered  and  converted  into  electrical  energy,  which  is 
stored  in  the  vehicle  battery  or  other  energy  storage  devices  for 
vehicle  electronics  equipment,  improving  the  fuel  economy  of 
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vehicles.  Studies  conducted  by  Stobart  and  Milner  [3]  have  shown 
that  based  on  the  assumption  of  no  increase  in  mass,  a  current 
market  price  of  approximately  $1000  would  see  the  investment 
break  even  over  the  lifetime  of  the  vehicle  when  a  TEG  device  is 
utilised. 

Although  TEG  has  broad  application  prospects,  due  to  the  rela¬ 
tively  low  thermal  conversion  efficiency,  it  is  still  only  at  the  stage 
of  theoretical  research  and  experimental  demonstration  without 
large-scale  practical  applications.  The  conversion  efficiency  of  TEG 
can  be  improved  from  the  following  two  aspects:  increasing  the 
thermoelectric  material  optimal  value  of  ZT,  as  well  as  increasing 
the  temperature  difference  between  the  hot  and  the  cold  sides  of 
thermoelectric  module  (TEM)  [4].  On  the  one  hand,  since  the 
optimal  value  of  ZT  is  an  inherent  property  of  thermoelectric  ma¬ 
terial,  improving  the  structure  of  thermoelectric  material  itself  as 
much  as  possible  can  raise  ZT  value.  But  ZT  values  of  large-scale 
commercial  thermoelectric  materials  are  generally  not  high,  for 
example,  this  value  of  B^Tes  is  only  about  1.  On  the  other  hand,  if 
conditions  permit,  efforts  can  be  made  on  the  optimization  of  TEG 
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heat  exchanger  structure,  such  as  improving  the  thermal  environ¬ 
ment  to  reduce  heat  loss,  making  the  contact  surface  as  flat  as 
possible  to  reduce  the  thermal  contact  resistance,  reducing  the 
impact  of  intermediate  links,  etc.  Increasing  the  temperature  dif¬ 
ference  across  TEMs  can  also  improve  the  conversion  efficiency  of 
TEG. 

Recently,  for  the  purpose  of  improving  the  output  power  and  the 
conversion  efficiency  of  TEG,  researchers  have  conducted  lots  of 
work  on  the  numerical  analysis  and  optimization  of  TEG  heat 
exchanger  models.  Esarte  et  al.  [5]  made  an  analysis  of  the  influence 
of  fluid  flow  rate,  heat  exchanger  geometry,  fluid  properties  and 
inlet  temperatures  on  the  power  supplied  by  TEG,  these  results  can 
give  those  who  have  to  design  TEG  a  good  idea  about  which 
operating  conditions  best  meet  the  specifications  required  for  a 
particular  application.  Suzuki  and  Tanaka  [6,7]  designed  two  types 
of  geometric  structure  for  TEG,  multi-panels  and  cylindrical  multi¬ 
tubes.  The  two  structure  types  were  both  exposed  to  two  thermal 
fluids  directly,  and  the  electric  power  generated  by  TEG  was  esti¬ 
mated  respectively.  The  output  power  of  the  proposed  15  systems 
for  multi-panels  were  deduced  from  heat  transfer  theory,  and  6 
systems  for  cylindrical  multi-tubes  were  also  investigated  in  that 
way  of  work.  The  results  showed  that  each  type  of  structure  had  the 
maximum  power  output  in  the  case  of  ideal  isothermal  systems, 
and  in  the  other  realistic  systems,  there  existed  a  certain  system 
that  supplied  the  maximum  power  output.  Meng  et  al.  [8]  studied 
the  complete  model  of  a  TEG  with  finned  heat  exchanger,  taking 
into  account  the  thermal  irreversible  effects.  Through  the  simula¬ 
tion,  it  showed  that  TEMs  can  obtain  0.13  W  output  power  and 
0.87%  the  maximum  conversion  efficiency,  it  had  also  pointed  out 
that  the  modelling  method  can  effectively  be  applied  to  the  low- 
quality  heat  recovery  TEGs.  Crane  and  Jackson  [9]  investigated 
thermoelectric  waste  heat  recovery  for  current  thermoelectric 
materials  with  advanced  heat  exchangers.  A  numerical  model  for 
heat  exchanger  integrated  with  thermoelectric  models  were 
created  and  validated  against  experimental  data.  It  showed  that  the 
heat  exchangers  can  achieve  power  densities  over  40  W/l.  Wang 
et  al.  [10]  presented  a  mathematical  model  of  a  TEG  device  using 
the  exhaust  gas  of  vehicles  as  heat  source.  The  model  simulated  the 
influence  of  the  relevant  factors  on  the  output  power  and  efficiency, 
such  as  exhaust  mass  flow  rate,  temperature  and  mass  flow  rate  of 
different  types  of  cooling  fluid,  convection  heat  transfer  coefficient, 
etc.  Hsiao  et  al.  [11]  established  a  numerical  model  combined  a 
TEM,  a  cooling  system  and  a  one-dimensional  thermal  resistance 
model,  and  applied  this  model  on  two  different  potential  positions, 
the  exhaust  tube  and  the  radiator,  to  examine  the  feasibility.  The 
maximum  power  density  produced  by  the  TEG  was  51.13  mW/cm2 
when  the  temperature  difference  is  290  °C.  The  results  showed  that 
TEM  presented  better  performance  on  the  exhaust  tube  than  on  the 
radiator. 

Currently,  TEGs  that  being  installed  on  the  vehicle  exhaust  tube 
mainly  have  the  structure  of  flat  panels.  Fig.  la  and  b  [12,13] 
respectively  show  the  single-layered  panel  and  the  multi-layered 
panel  structure.  The  panel-structured  heat  exchangers  have  the 
main  advantages  of  simple,  stable  and  reliable,  making  the  TEMs 
heated  and  cooled  sufficiently,  lower  cost  and  so  on.  However,  the 
panel-structured  heat  exchanger  changes  the  shape  of  the  exhaust 
tube,  which  has  lower  portion  of  the  space  under  vehicle  chassis, 
causing  difficulties  in  installation  and  maintenance.  When  the  ge¬ 
ometry  shape  of  the  exhaust  tube  changes,  exhaust  back  pressure 
will  be  affected  to  some  extent,  which  leads  to  the  influence  of  the 
engine  operating  conditions.  These  factors  limit  large  scale  appli¬ 
cation  of  the  panel  structured  heat  exchanger  TEG. 

Based  on  summarizing  the  advantages  and  disadvantages  of 
conventional  structures  of  TEG,  this  paper  proposes  a  newly 
designed  TEG  with  cylindrical  shell  and  straight  fins.  Since  the  heat 


Fig.  1.  a.  TEG  with  single-layered  panel  heat  exchanger,  b.  TEG  with  multi-layered 
panel  heat  exchanger. 


transfer  process  during  the  operation  of  TEG  is  complex,  including 
the  heat  conduction  among  various  components  and  the  heat 
convection  between  fluids  and  solids,  the  heat  transfer  numerical 
models  for  heat  exchanger  and  TEM  are  established  separately, 
which  can  evaluate  the  performance  of  TEG  effectively.  The  goals  of 
the  numerical  simulations  include  the  temperature  changes  of 
exhaust  and  coolant  in  the  axial  direction,  the  temperature  distri¬ 
bution  on  the  fins,  and  the  variations  of  output  power  and  con¬ 
version  efficiency  of  TEG.  The  simulations  also  focus  on  the 
thickness  of  the  insulation  layer  of  TEG,  the  exhaust  and  the  coolant 
inlet  temperatures  as  well  as  the  inlet  mass  flow  rate,  the  TEM 
positioning  and  other  external  factors,  which  may  affect  the  TEG 
output  characteristics.  All  these  simulations  are  completed  under 
steady-state  conditions.  The  numerical  model  simulation  analyses 
can  help  the  TEG  designers  making  choices  of  the  relevant  pa¬ 
rameters,  and  provide  appropriate  guidance  to  complete  the 
designing  optimization  of  TEG.  It  also  lay  a  solid  foundation  for 
future  TEG  simulations  under  transient  conditions  of  vehicle 
engine. 

2.  Descriptions  on  the  structure  of  TEG  with  cylindrical  shell 
and  straight  fins 

The  conventional  TEG  which  is  on  vehicle  exhaust  tube  causes 
significant  structural  changes,  insufficient  contact  of  power  gen¬ 
eration  devices,  obvious  influence  on  the  engine  exhaust  back 
pressure,  etc.  To  overcome  these  defects,  this  paper  proposes  a 
newly  designed  TEG  with  cylindrical  shell  and  straight  fins,  whose 
schematic  is  shown  in  Fig.  2.  This  device  is  compact  in  structure  and 
has  high  operational  reliability,  it  can  be  also  arranged  on  the  po¬ 
sition  between  the  vehicle  exhaust  catalytic  converter  and  the 
silencer  to  make  the  effective  recovery  of  waste  heat.  Two  semi- 
cylindrical  shells  directly  contact  the  exhaust  tube,  which  are 
joined  as  one  unit  by  connecting  bolts.  Riveting  fasten  devices  are 
used  to  support  the  heat  transfer  fins,  TEMs  and  cooling  tubes.  The 
TEMs  and  heat  transfer  fins  are  in  direct  contact,  avoiding  the 
exhaust  tube  structural  transformation,  and  it  will  not  cause  any 
influence  on  the  engine  exhaust  back  pressure.  A  certain  number  of 
symmetry  fitting  grooves  on  the  semi-cylindrical  shells  are 
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supporting  disc  device 


Fig.  2.  Schematic  of  TEG  with  cylindrical  shell  and  straight  fins. 

reserved,  so  the  number  of  heat  transfer  fins  can  be  changed 
optionally,  accommodating  the  different  power  generation  needs. 
The  cooling  tubes  are  branches  of  the  engine  cooling  system,  and 
the  engine  coolant  flows  into  the  tubes  to  cool  down  the  cold  sides 
of  TEM.  The  cooling  tubes  are  flat  rectangular  cross-section  chan¬ 
nels,  which  have  unified  specifications  and  convenient  interfaces. 
Supporting  discs  with  bolts  connecting  the  semi-cylindrical  shells 
are  set  on  both  ends  of  the  cooling  tubes,  so  that  the  whole  TEG 
structure  can  be  hold  stably.  This  device  has  high  thermal  conver¬ 
sion  efficiency,  relatively  low  thermal  resistance,  and  it  can  improve 
the  exhaust  waste  heat  recovery  efficiency  compared  to  those 
conventional  vehicle  exhaust  TEG  systems.  It  is  also  suitable  for 
mass  production  and  usage. 

The  semi-cylindrical  shells  that  directly  contact  the  exhaust 
tube  are  made  of  stainless  steel,  which  is  able  to  meet  the  needs  of 
exhaust  temperature  for  heat  exchangers  (the  temperature  of 
exhaust  that  is  near  the  vehicle  catalytic  converter  is  about  500  °C). 
The  heat  transfer  fins  are  made  of  aluminium  alloy,  which  has  good 
thermal  performance  and  can  maintain  high  enough  temperature 
on  the  hot  sides  of  TEM.  A  certain  amount  of  thermal  grease  are 
smeared  between  the  heat  transfer  fins  and  the  hot  sides  of  TEM, 
reducing  the  thermal  contact  resistance  caused  by  the  surface 
roughness  in  order  to  improve  the  thermal  conductivity.  The 
cooling  tubes  are  also  made  of  aluminium  alloy.  Due  to  the  straight 
channel  and  thin-walled  rectangular  cross-section  structure,  they 
can  directly  be  welded  by  using  aluminium  alloy  sheets.  The 
aluminium  alloy  cooling  tubes  are  fixed  on  the  TEG  through  two 
supporting  discs  with  fastening  devices  at  both  ends.  A  certain 
preload  is  made,  so  the  TEMs  and  the  heat  transfer  fins  can  contact 
closely.  Between  the  exhaust  tube  and  the  supporting  discs,  as  well 
as  the  supporting  discs  and  the  connecting  bolts  of  the  semi- 
cylindrical  shells,  there  exists  some  special  washers,  which  are 
made  of  insulation  material  and  can  minimize  the  influence  of  the 
supporting  discs  heating  effect  on  the  coolant. 

3.  Heat  transfer  numerical  model  and  validation 

3.1.  Numerical  model  description 

From  the  schematic  of  TEG  with  cylindrical  shell  and  straight 
fins  which  is  shown  in  Fig.  2,  it  can  be  seen  that  the  exhaust  flows 
through  the  tube,  so  there  exists  heat  convection  between  the 
exhaust  and  the  tube  wall.  After  being  conducted  through  the  tube 
wall  and  the  cylindrical  shell  axially  and  radially,  one  part  of  heat 
flux  flows  to  the  heat  transfer  fins  and  the  TEMs,  while  the  other 
part  is  released  into  the  surrounding  environment  by  natural 
convection.  The  heat  flux,  which  is  conducted  to  the  TEMs,  is 
divided  into  two  parts:  one  is  transferred  as  electric  power,  and  the 
other  is  conducted  to  the  cold  sides  of  TEM. 

According  to  the  general  theory  of  numerical  heat  transfer,  the 
first  work  to  be  carried  out  in  modelling  is  the  regional 


discretization  [14].  As  it  is  shown  in  Fig.  3  (due  to  the  symmetrical 
distribution  of  the  heat  transfer  fins,  there  is  only  one  pair  of  fins  in 
the  figure  for  simplification,  the  insulation  layer  and  the  cylindrical 
shell  are  also  simplified  to  one  layer),  segmentation  is  made  on  the 
cylindrical  shell,  the  TEMs  and  the  cooling  tubes  in  axial  direction. 
The  control  volumes  (CVs)  are  uniform  in  length,  which  is  Ax 
respectively,  and  i  represents  the  number  of  CV.  Since  TEG  is  central 
symmetry  about  the  axis,  conditions  are  all  the  same  along  the 
circumferential  direction.  The  internal  node  method  is  used  in  grid 
generating  process.  The  exhaust  flow,  exhaust  tube  wall,  cylindrical 
shell,  cooling  tube  wall  and  coolant  are  all  in  one-dimensional  axial 
grid  structure.  To  research  the  TEM  output  performance  on 
different  locations  of  the  heat  transfer  fins,  the  thermocouples  and 
the  heat  transfer  fins  involved  in  segment  i  along  the  axial  direction 
are  divided  into  radial  direction  grids,  whose  lengths  are  Ay,  so  a 
two-dimensional  grid  structure  is  established  on  the  heat  transfer 
fins.  The  following  assumptions  are  made  to  simplify  the  complex 
problem  of  modelling  [1]: 

i.  The  axial  heat  conduction  within  the  thermocouples  is 
ignored  due  to  the  dominion  of  the  transverse  conduction 
along  the  thermocouples. 

ii.  Thermal  resistance  through  the  exhaust  tube  wall,  the  cy¬ 
lindrical  shell,  the  insulation  layer,  the  ceramic  plates,  the 
metallic  strips  of  TEMs  and  the  cooling  tube  wall  is  taken  into 
account  as  an  equivalent  thermal  contact  resistance. 

iii.  The  radial  length  of  the  two-dimensional  CV,  Ay,  is  larger 
than  the  length  of  one  single  thermocouple,  so  the  integrity 
of  thermocouple  can  be  ensured  to  complete  the  analysis. 

iv.  The  gap  between  thermocouples  is  ignored. 

v.  Ignoring  the  radiation  heat  transfer  and  the  heat  losses  be¬ 
tween  thermocouples  and  the  heat  transfer  fins  or  the 
cooling  tubes. 

vi.  The  conducting  heat  flux  between  the  insulation  layer  and 
the  heat  transfer  fins  is  ignored. 

Based  on  these  assumptions,  heat  balance  equations  can  be 
obtained  for  each  CV  of  TEG  numerical  model. 

The  heat  balance  equation  for  exhaust  flow  is: 

Qcon_ex,i  =  ^ex^ex  (Jex,i  ~  ^exj+l ) 

—  ^ex-wall^^ex_tube^(^ex,i  —  ^tw,i)  O) 

The  heat  transfer  coefficient  between  the  exhaust  flow  and  the 
exhaust  tube  wall  hex-waii  is  given  by  the  following  equations  1]: 

^ex-wall  =  NtJex^ex/tfex_tube 


(2) 
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(ff/2)  (Ref  -  lOOo)  Prf 
l+2.7(/f/2)1/2(Pr^3-l) 


(3) 


The  heat  balance  equations  of  the  insulation  layer  in  each  CV 
are: 

Qcyl, radial  **  Qins, radial  +  Qins, axial  +  ^fin^finO  (7-1) 


ft 


A  + 


B 


(4) 


Qins,  radial  — 


^ins  ,i  ^surf  ,i 


1 


27rAinsAx 


ln(l  + 


_ 2  tjns _  j 

dtube  +2  ttube  +2  tcy|  ^ 


(7.2) 


7rAin 


Qins, axial  — 


(  2be  +  ftube  +  fCyl  +  fins)  (  2be  +  fftibe  +  fcyl) 


(2^ins,i 


ins,i+l  ^ins,i-l 


Ax 


(7.3) 


where  for  laminar  flow  {Ref  <  2100),  A  =  0,  B  =  16.0  and  m  =  1.0,  for 
transition  flow  (2100  <  Re f  <  4000),  A  =  0.0054,  B  =  2.3  x  10  8  and 
m  =  -0.6667,  and  for  turbulent  flow  {Re f  >  4000),  =  0.00128, 

B  =  0.1143  and  m  =  3.215. 

The  exhaust  tube  wall,  the  cylindrical  shell  and  the  heat  insu¬ 
lating  layer  can  be  viewed  as  three  cylindrical  layers  (ignoring  the 
volume  of  the  insulation  layer  being  spaced  among  the  heat 
transfer  fins),  heat  conduction  process  exists  along  the  axial  and  the 
radial  directions.  The  heat  flux  which  is  conducted  from  the  inter¬ 
nal  layer  along  the  radial  direction  can  be  thought  as  the  total  input 
heat  flux  to  the  external  layer,  so  the  following  four  sets  of  heat 
balance  equations  can  be  obtained. 

The  heat  balance  equations  of  the  exhaust  tube  wall  in  each  CV 
are: 

Qcon  _ex,i  =  Qex_tube,  radial  +  Qex_tube,  axial  (5-1) 


Qex_tube,  radial 


^tw.i  Tcyl.i 

_J _ In  ( 1  +  2fex_tube\ 

ex_tube^^  \  ^ex_tube  J 


(5.2) 


qfin0  is  the  heat  flux  being  conducted  to  one  of  the  heat  transfer  fins 
in  each  CV,  Nfin  is  the  total  number  of  the  heat  transfer  fins. 

The  heat  balance  equations  between  the  surface  of  insulation 
layer  and  the  external  environment  in  each  CV  are: 


Qconjns  —  ^con_air^(  ^ex_tube  +  2  tex  tube  +  2fCyi  +  2tjr 


Ax(rs 


surf,i 


-  T  ■ 

1 


(8.1) 


Qconjns  —  Qins, radial  (8.2) 

The  natural  air  convection  coefficient  on  the  surface  of  the 
insulating  layer  hcon_ air  is  [15]: 

^con_air  =  ^wcon_air^air  /  (  ^ex_tube  +  ^fex_tube  +  2 tCyi  +  2tjns  J  (9) 


Nucon_air  =  C(Grair  Prair)n  =  C{Ra)n 


(10) 


TZX. 


ex_tube 


Qex_tube,  axial  — 


(dex_tube  ,  f.  \  ^  /  deX_tube  \  Z 

2  rex_tube  J  ~  ^  2  J 


(2Ttw,i 


tw,i+l  ^tw,i-l 


Ax 


(5.3) 


The  heat  balance  equations  of  the  cylindrical  shell  in 


are: 


Qex_tube,  radial  —  Qcyl,  radial  +  Qcyl,  axial 
^cyl,i  —  ^ins  ,i 


For  the  horizontal  cylinder,  when  10 4  <Ra<  109,  C  =  0.53,  n  =  1/4, 
each  CV  when  ^  <Ra<  1()12’  C  =  ai3>  n  =  V3- 

7<§t^ins  —  ^air)  f^tube  +  ^ftube  +  2 tCyl  +  2 fins) 

(6.1)  Crair  = - ^ - 2 - - (11) 


Qcyl,  radial  —  ^ 


In  1  + 


2  tcv 


2TrAcyi  Ax  y  ^eX  _tube  +2  fex_tube 


(6.2)  2 

7  = 


T  +  T 

1  ins  ^  1  air 


(12) 


tuA( 


cyl 


Qcyl,  axial  — 


(%^  +  texjtube  +  tCyl)'  -  +  Wtube)2]  p :yl,i  -  Tc yU+1  -  Tc yW,) 


Ax 


(6.3) 
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For  the  heat  transfer  fins,  the  schematic  of  discretization  along 
the  radial  direction  is  shown  in  Fig.  4.  Where,  b  is  the  radial  distance 
from  the  root  of  the  heat  transfer  fin  to  TEM,  a  is  the  side  length  of 
TEM. 

Since  the  heat  transfer  fins  are  equipped  with  TEMs,  the  heat 
transfer  conditions  are  not  all  the  same  for  different  positions.  The 
thermal  analyses  of  the  fins  should  be  divided  into  three  parts.  For 
0  <y<b,  there  exists  natural  air  convection  on  the  surface  of  the 
heat  transfer  fin,  the  temperature  distribution  is  in  the  presence  of 
the  second-order  ordinary  differential  equation: 


d2rfln 

dy 2 


2frfin_air 

^fin^fin 


(13.1) 


The  central  difference  of  second-order  accuracy  is  used  to 
transfer  the  ordinary  differential  equation  into  the  following  form: 


^fin,i(/<+l)  2Tfmik  +  7fin,i(/<-l) 

Ay2 


2ftfin_air,t7< 

^fin^fin 


(13.2) 


The  initial  conditions  are: 


R  —  2aLj£ieg  j  AfEieg  +  ^copper  (18) 

where  LTmeg  and  ATEieg  are  the  length  and  the  cross  sectional  area  of 
the  thermocouple  leg  respectively. 

All  the  thermocouples  are  connected  in  series,  so  the  current  /  is 
determined  by  the  number  of  thermocouples: 


/  m 


n  k2 

S  a  (^th,ik  — 
i=  1  k  =  \ 


^tc,i/c)^fin 


RNjE  +  Rl 


(19) 


Ri  is  the  total  external  resistance.  The  heat  absorbed  by  the 
thermocouple  in  the  whole  CV  is: 


k2 

*7heat_th,i  =  Qh,ik  (20) 

k  =  1 

For  a  +  b  <  y  <  Hfin,  the  temperature  distribution  equation  is  the 
same  with  Eqs.  (13.1)  and  (13.2),  so  the  total  quantity  of  heat 
convention  between  the  fins  and  the  air  is: 


^ns,4nAxttin7flnl2Ay7finn 

<7fin0  -  2/lfin_air,ilAx4y(7fin,ii  -  Cair) 


(14) 


9air_convection2 


2  J2  htm_^AlM^y(TtmJk  ~  fair) 


(21) 


where  the  hfm_air,ik  is  the  ikth  CV’s  two-dimensional  natural  con¬ 
vection  coefficient.  The  total  quantity  of  heat  convection  on  one  fin 
is  (the  upper  and  the  lower  side  of  the  fin  shown  in  Fig.  4  both  has 
the  same  quantity  of  natural  convention,  and  ignoring  the  area  of 
two  vertical  sides  since  it  is  tiny): 


In  Eqs.  (15)  and  (19)-(21),  k\,  k2,  k3  refer  to  the  different 
numbers  of  discrete  units  in  each  interval  of  y. 

At  last,  the  heat  balance  equation  of  a  single  fin  in  axial  CV  is: 

*7fin0  ~  ^air_convectionl  “F  9heat_th,i  "F  <7air_convection2  (22) 


*7air_convectionl  —  2  hfin_alr, ik^xhy  ( Tfm  jk  ^airj  O^) 

k  =  1 

For  b  <  y  <  a  +  b,  there  exists  the  similar  temperature  distri¬ 
bution  second  ordinary  differential  equations: 


d27fin  =  Qh^TE 
dy2  ^fin^fin^^fin 


(16.1) 


^fin,i(/<+l)  2Tfini/<  +  1 )  _  Qh,ik^TE 

Ay2  ^fin^fin^a^fin 

In  Eq.  (16.2): 

Qh.ik  =  ^ th.il J  +  ,ik  ~  ^tc ,ik)  2!  ^ 


(16.2) 


(17) 


where  Tth,ik  is  equal  to  Tfin, iio  R  is  the  internal  resistance  of  one 
thermocouple,  which  is  the  sum  of  the  resistance  of  the  thermo¬ 
couple  legs  and  the  copper  strip,  that  is: 


In  the  CV  shown  in  Fig.  3,  the  heat  released  from  the  cold  side  of 
TEM  is: 


Qheat_tc,i  —  Qc,ik 
k  =  0 


(23.1) 


Qc,ik  = 


^tc.ilJ  “F  Tfc  ik)  +  2^  R 


NjE^xAy 

Nf[nLci 


(23.2) 


The  heat  balance  equations  of  the  cooling  tube  wall  in  each  CV 


are: 


Qheat_tc  —  Qcw, radial  “F  Q cw, axial 

(24.1) 

^  ]  Avilf  Ttc,i  —  Tcti 

Qcw,  radial  —  ^ctAXWct  ^ 

(24.2) 

^  ^ctWcttct(2 Tcti  -  Tct  j+1  -  Tct,i-i) 

9 cw,  axial  —  ^ 

(24.3) 

Ttc, i  is  the  mean  value  of  Ttc,ik  on  the  whole  discretization  zone. 
The  heat  balance  equation  between  the  coolant  and  the  cooling 

tube  wall  is: 

*7con_cf,i  —  ^cf  Qf  (TcU  ^cf,i+l)  —  ^cf-wallWctAx^ 

=  9cw, radial 


Trti~Tf 


cf  ,i 


The  output  power  of  TEG  is: 


(25) 
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P  —  f^heat_th.i  9heat_tc,i')^fin  (26) 

i'  =  l  i'  =  l 

The  TEM  conversion  efficiency  is: 

^TEM  =  ~h -  (27) 

X]  9heat_th^fm 
i=l 

and  the  TEG  conversion  efficiency  of  the  whole  numerical  model  is: 

^TEG  =  ~n -  (28) 

X]  Qcon_  ex,i 
i  =  l 


3.2.  Numerical  process  and  validation 

Along  the  axial  direction,  Ax,  which  is  the  length  of  each  CV, 
ranges  from  20  to  40  cm.  Due  to  the  mild  temperature  gradients 
across  the  exhaust  tube,  this  level  of  discretization  is  adequate 
enough  for  the  accurately  modelling  of  the  TEG  numerical  model 
across  the  full  exhaust  tube  length.  Because  the  radial  dimensions 
of  the  fins  are  relatively  small,  the  radial  direction  length  of  each  CV 
Ay  varies  from  0.4  to  1  cm.  The  solution  methodology  of  this  nu¬ 
merical  model  is  implemented  by  using  an  iterative  Newton- 
Raphson  method.  Since  all  the  TEMs  on  one  fin  is  connected  elec¬ 
trically  in  series,  the  current  I  is  constant  for  all  axial  positions, 
which  may  influence  the  iteration  speed  of  Newton-Raphson 
method,  but  it  is  still  in  a  reasonable  range. 

In  order  to  validate  the  accuracy  of  this  numerical  model,  a 
finite  elements  model  is  built  and  analysed  in  ANSYS  Workbench 
environment.  The  number  of  fins  is  8,  and  all  of  the  fins  are 
symmetrical  along  the  radial  direction,  so  it  is  more  practical  and 
resource-saving  to  take  one  fin  to  simulate.  One-eighth  of  the 
structure  of  this  model  is  shown  in  Fig.  5.  The  boundary  of  the 
exhaust  gas  domain  is  added  with  the  inflation  layer  to  adapt  the 
complex  conditions  on  the  interface  between  solid  and  fluid. 
There  are  5  layers  of  the  inflation  and  the  total  thickness  is  1  cm. 
Other  boundary  conditions  of  this  finite  elements  model,  such  as 
the  natural  convention  between  the  outside  surfaces  of  TEG 
components  and  air,  the  inlet  and  the  outlet  parameters  of  the 
exhaust  and  the  cooling  fluid,  are  consistent  with  the  numerical 
model. 

Figs.  6  and  7  show  the  output  power  and  the  heat  conversion 
efficiency  comparisons  between  the  two  models  respectively  un¬ 
der  certain  conditions  when  the  exhaust  inlet  temperature 
changes.  It  can  be  observed  that  there  are  tiny  discrepancies  be¬ 
tween  the  two  models,  mainly  due  to  the  intrinsic  differences  of 
the  implementation  of  the  two  programs.  The  assumptions  i  and  iv 
mentioned  above  could  not  be  supported  in  ANSYS  simulations, 
because  the  axial  heat  conduction  of  TEM  and  the  gaps  between 
them  are  taken  into  consideration  in  the  program,  which  may 


Fig.  5.  Schematic  of  the  finite  elements  model  structure  with  grids  (one-eighth). 


Fig.  6.  Comparisons  of  the  output  power  between  the  numerical  model  in  Matlab  and 
the  finite  elements  model  in  ANSYS. 


influence  the  output  performance  of  the  ANSYS  model.  Besides, 
with  the  increasing  number  of  CV,  the  calculation  precision  of  the 
numerical  model  can  be  improved  as  well,  leading  to  better 
consistence  between  these  two  models.  Despite  of  these,  the  dif¬ 
ferences  between  the  two  models  are  within  the  allowable  range, 
which  validates  that  the  proposed  numerical  model  could  be  an 
equivalent  thermal  simulator  to  ANSYS  in  thermoelectric 
simulation. 


4.  Simulation  results  and  discussion 

In  the  following  numerical  simulations,  only  under  the  steady 
state  working  conditions  the  TEG  output  performances  are  ana¬ 
lysed.  Transient  working  conditions  analyses,  which  can  study  the 
output  performances  of  TEG  under  the  whole  working  cycle  of 
engine,  should  be  combined  with  a  special  engine  simulation 
software,  such  as  GT-POWER,  for  example,  and  this  is  not  in  the 
scope  of  this  discussion. 


Fig.  7.  Comparisons  of  the  TEM  and  the  TEG  conversion  efficiencies  between  the  nu¬ 
merical  model  in  Matlab  and  the  finite  elements  model  in  ANSYS. 
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4  A.  Parameter  settings 

Choosing  a  model  with  diesel  engine  exhaust  tube  size  that 
matches  the  numerical  simulation.  The  inner  diameter  of  the 
exhaust  tube  is  60  mm,  and  the  wall  thickness  is  3  mm.  The  length 
of  the  cylindrical  shell  is  685  mm,  and  the  thickness  is  10  mm.  In 
order  to  ensure  the  heat  losses  as  little  as  possible,  and  raise  the 
temperature  of  the  hot  side  of  TEG,  there  sets  an  insulation  layer 
outside  the  cylindrical  shell.  The  radial  height  of  the  heat  transfer 
fin  is  60  mm,  the  thickness  is  5  mm,  and  the  length  is  equal  to  that 
of  the  cylindrical  shell.  The  exhaust  tube  and  the  cylindrical  shell 
are  made  of  cast  iron  and  nickel-chromium  steel  respectively,  and 
the  material  of  the  insulation  layer  and  the  heat  transfer  fin  are 
refractory  slag  wool  and  aluminium-silicon  alloy  respectively.  The 
cooling  tube  is  also  made  of  aluminium-silicon  alloy  material,  the 
cross-sectional  dimensions  are  34  mm  x  10  mm,  and  the  thickness 
is  1  mm. 

As  the  engine  exhaust  temperature  is  relatively  high,  the 
semiconductor  thermoelectric  material  in  the  middle  temperature 
region  (450-850  K)  PbTe  is  used  as  the  material  of  TEM,  the 
dimension  of  each  thermocouple  is  1.5  mm  x  1.5  mm  x  2.0  mm, 
and  the  relevant  physical  parameters  are  a  =  2.4  x  10  4  V/K, 
A  =  0.5  W/(K  m)  ,  <7  =  1.2  x  10-5  Q  m  [16].  The  P  and  N  junctions  of 
the  thermocouples  are  connected  by  copper  strips  with  the  thick¬ 
ness  of  0.3  mm.  Between  the  copper  strips  and  the  heat  transfer 
fins,  as  well  as  the  copper  strips  and  the  cooling  tubes,  there  exists 
some  ceramic  plates  with  the  thicknesses  of  0.7  mm.  Although  the 
thermal  physical  parameters  of  the  material  change  with  the 
temperature,  in  order  to  facilitate  the  numerical  simulation,  the 
temperature  which  is  near  the  base  of  the  heat  transfer  fins  is 
selected  as  the  reference  temperature.  In  the  entire  simulation 
process,  the  physical  parameters  are  set  to  be  constant.  The  load 
resistance  value  is  set  to  be  equal  to  the  total  internal  resistance  of 
the  thermocouples,  so  the  maximum  output  power  can  be  achieved 
under  the  same  conditions  [17].  Tables  la  and  b  show  the  relevant 
physical  parameters  in  numerical  simulation. 

The  inlet  temperature  of  the  exhaust  flow  is  set  to  be  773  K 
(500  °C),  and  the  flow  velocity  is  60  m/s.  As  the  coolant  in  the  TEG 
cooling  system  is  supplied  by  the  engine  cooling  system,  the  inlet 
temperature  of  the  coolant  is  set  to  be  353  I<  (80  °C),  and  the  flow 
velocity  is  1  m/s.  In  different  simulation  conditions,  the  inlet  tem¬ 
perature  and  the  flow  velocity  of  the  exhaust  and  the  coolant  are 
recognized  as  the  adjustable  variables.  The  correlations  mentioned 
above  on  the  heat  transfer  coefficient  in  the  numerical  model  can 
still  be  chosen  because  the  Reynolds  and  Prandtl  numbers  are  still 
within  the  extent  of  its  validation.  Exhaust  and  coolant  flows  are 
parallel  flow  type. 

4.2.  Simulation  results 

Fig.  8  shows  the  temperature  variations  of  the  exhaust  and  the 
coolant  along  the  axial  direction  when  the  TEG  model  is  axially 
divided  into  20  average  discrete  CVs,  the  thickness  of  the  insulation 


Table  la 

Relevant  physical  parameters  of  solids  in  the  numerical  simulation. 


P  (kg/m3) 

X  (W/(m  I<)) 

a  m) 

a  (V/K) 

Exhaust  tube 

7570 

28.7 

- 

- 

Cylindrical  shell 

7820 

22.2 

- 

- 

Insulation  layer 

350 

0.067 

- 

- 

Heat  transfer  fin 

2660 

180 

- 

- 

Thermocouple 

8200 

0.5 

1.2  x  10"5 

2.4  x  10"4 

Copper  strip 

8930 

372 

1.75  x  10"8 

- 

Ceramic  plate 

3970 

17 

1.1  x  102 

- 

Cooling  tube 

2660 

180 

- 

- 

Table  lb 

Relevant  physical  parameters  of  liquid  and  gases  in  the  numerical  simulation. 


P  (kg/m3) 

X  (W/(m  K)) 

CP  (J/(kg  K)) 

P  (kg/(m  s)) 

Pr 

Exhaust  flow 

0.457 

0.0656 

1185 

3.48  x  10"5 

0.63 

(773  K) 

Coolant  (353  I<) 

971.8 

0.674 

4195 

3.551  x  10"4 

2.21 

Air  (300  K) 

1.165 

0.0276 

1005 

1.86  x  10"5 

0.701 

layer  is  10  mm  and  b  is  15  mm.  As  it  can  be  seen,  the  temperature  of 
exhaust  and  the  coolant  are  linear  variations.  Different  from  the 
conventional  heat  exchangers,  due  to  the  Seebeck  effect,  the  TEM 
transforms  part  of  the  heat  absorbed  from  the  hot  sides  into  elec¬ 
trical  energy,  resulting  the  heat  released  to  the  coolant  is  less  than 
that  absorbed  from  the  hot  sides  of  TEM.  So  the  temperature  dis¬ 
tribution  is  linear,  and  this  conclusion  is  consistent  with  the  liter¬ 
ature  [1].  Unlike  [1],  which  used  water  as  both  the  hot  and  the 
cooling  fluids  to  cause  the  similar  temperature  variation  ranges, 
during  the  simulation  process  in  this  paper,  the  hot  fluid  is  exhaust 
whose  specific  heat  capacity  and  density  are  relatively  small,  while 
the  cooling  fluid  is  engine  coolant  whose  specific  heat  capacity  and 
density  are  relatively  high,  so  the  temperature  variation  range  of 
the  exhaust  is  larger  than  that  of  the  coolant.  Due  to  this  limited 
coolant  temperature  variation  range,  there  is  not  too  much  influ¬ 
ence  on  the  performance  of  TEG  whether  the  exhaust  and  the 
coolant  flow  type  is  parallel  or  counter,  and  this  will  bring  great 
convenience  on  practical  applications  of  engineering. 

Fig.  9  shows  the  temperature  distributions  along  the  radial  di¬ 
rection  for  each  CV  on  one  fin  (all  the  fins  are  radially  symmetrical, 
so  it  takes  one  fin  to  represent  all  without  considering  the  influence 
of  convention  caused  by  the  gravity  of  air),  the  thickness  of  the 
insulation  layer  is  10  mm  and  b  is  15  mm,  the  natural  convection 
heat  transfer  coefficient  is  assumed  to  be  constant.  As  it  can  be  seen, 
the  temperature  distributions  of  the  fin  are  not  linear,  larger  tem¬ 
perature  drops  happens  near  the  base  of  the  fin.  In  the  middle  part 
of  the  fin,  which  is  also  the  place  of  the  TEM  hot  sides  being  laid, 
heat  flux  is  absorbed  by  TEM,  so  the  temperature  continues  to  drop. 
There  is  relatively  little  amount  of  heat  being  conducted  to  the  last 
part  of  the  fin,  and  it  is  all  taken  away  to  the  environment  through 
natural  convection,  the  temperature  drop  is  relatively  slow.  It  is 
concluded  that,  since  the  severe  temperature  drop  near  the  base  of 
the  fins,  the  arrangements  of  TEM  should  be  as  near  to  the  base  of 
the  fins  as  possible.  In  other  words,  the  smaller  b  is,  the  higher 
temperature  of  the  hot  sides  of  TEM  is.  The  related  discussions  will 
be  given  later  in  this  section. 


2  750 


353.2  2 


x(m) 


Fig.  8.  Temperature  variations  of  the  exhaust  and  the  coolant  along  the  axial  direction. 
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Fig.  10  shows  the  output  power,  the  TEM  conversion  efficiency 
and  the  TEG  conversion  efficiency  of  each  CV,  the  thickness  of  the 
insulation  layer  is  10  mm  and  b  is  15  mm.  As  it  can  be  seen,  with  the 
axial  temperature  drop  of  the  exhaust,  the  output  power  and  the 
conversion  efficiencies  of  each  CV  decrease  linearly.  The  TEG  con¬ 
version  efficiency  is  less  than  the  TEM  conversion  efficiency  about 
2%,  which  is  caused  by  the  heat  losses  through  natural  convection 
on  the  surface  of  the  insulation  layer  and  the  heat  transfer  fins. 

The  following  numerical  analysis  mainly  focus  on  the  adjustable 
variables  of  TEG.  To  evaluate  the  influence  degree  of  each  variable 
to  the  output  performance  of  TEG,  an  evaluation  index  is  defined  as 
“the  average  variable  contribution  rate  to  the  output  power”,  CRvan 
which  is  in  the  following  form: 


CRvar  — 


_ Pm+ 1  -Pm _ 

^  variablem+i  -variable™ 
m  =  1 


j- 1 


(29) 


where  j  refers  to  the  number  of  the  sample  points  of  the  variable 
variation  range,  Pm  and  Pm+i  refer  to  the  output  power  of  the  whole 
TEG  model  at  the  adjacent  sample  points,  meanwhile  variablem  and 
variablem+i  refer  to  the  corresponding  value  of  variable  at  the 
adjacent  sample  points. 


Fig.  10.  Output  power,  TEM  conversion  efficiency  and  TEG  conversion  efficiency  of 
each  CV. 


Firstly,  assuming  the  material  of  insulation  layer  is  unchanged, 
the  output  power,  the  TEM  conversion  efficiency  and  the  TEG 
conversion  efficiency  of  each  CV  are  researched  under  the  condi¬ 
tions  that  the  thickness  of  insulation  layer  is  5  mm,  10  mm  and 
15  mm  respectively,  b  is  15  mm,  as  showed  in  Fig.  11a— c.  It  is  seen 


Fig.  11.  a.  Output  power  of  each  CV  under  different  thicknesses  of  the  insulation  layer, 
b.  TEM  conversion  efficiency  of  each  CV  under  different  thicknesses  of  the  insulation 
layer,  c.  TEG  conversion  efficiency  of  each  CV  under  different  thicknesses  of  the 
insulation  layer. 
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from  the  figures  that  the  output  performances  of  each  CV  drops 
linearly,  which  is  consistent  with  the  conclusion  obtained  from 
Fig.  10.  The  thicker  the  insulation  layer  is,  the  higher  output  power 
and  conversion  efficiencies  are.  It  is  because  that  the  thermal 
resistance  increases  with  the  thickness  of  the  insulation  layer,  so 
the  heat  conducted  through  the  insulation  layer  becomes  less,  and 
the  heat  being  taken  away  to  the  environment  through  air  natural 
convection  becomes  less  too.  The  total  amount  of  the  heat 
exchanged  between  the  exhaust  flow  and  the  tube  wall  is  the 
summation  of  the  heat  conducted  to  heat  transfer  fins  and  those 
being  taken  away  through  the  air  natural  convection  on  the  surface 
of  the  insulation  layer,  so  increasing  the  thickness  of  the  insulation 
layer  helps  to  enhance  heat  conducting  to  the  heat  transfer  fins,  the 
output  power  and  the  conversion  efficiencies  can  naturally  be 
higher  as  well.  The  output  power  for  the  whole  TEG  model  should 
be  recognized  as  the  Pm  in  Eq.  (29),  the  average  insulation  layer 
thickness  contribution  rate  to  the  output  power  CRinsuiation  layer 
thickness  is  0.49  W/mm. 

Fig.  12  shows  the  variations  of  output  power,  the  TEM  conversion 
efficiency  and  the  TEG  conversion  efficiency  for  the  whole  TEG 
model  with  the  inlet  temperature  of  exhaust,  the  thickness  of  the 
insulation  layer  is  10  mm  and  b  is  15  mm.  As  it  can  be  seen,  with  the 
increasing  of  the  exhaust  inlet  temperature,  more  heat  is  carried  to 
conduct  through  the  exhaust  tube  wall,  so  the  output  performances 
improve  as  well.  The  average  exhaust  inlet  temperature  contribution 
rate  to  the  output  power  CRexhaust  inlet  temperature  is  0.52  W/K.  In  the 
engineering  applications,  TEG  can  be  arranged  on  the  front  position 
of  the  exhaust  tube  as  possible  to  get  a  higher  exhaust  temperature. 

In  addition  to  improve  the  exhaust  inlet  temperature,  it  can  also 
reduce  the  coolant  temperature  to  enhance  the  output  perfor¬ 
mance  of  TEG,  both  of  these  methods  mean  to  increase  ATTe-  Fig.  13 
shows  the  variations  of  the  output  power,  the  TEM  conversion  ef¬ 
ficiency  and  the  TEG  conversion  efficiency  for  the  whole  TEG  model 
with  the  inlet  temperature  of  the  coolant,  the  thickness  of  the 
insulation  layer  is  10  mm  and  b  is  15  mm.  It  is  clear  that  as  the 
coolant  temperature  decreases,  the  output  power  and  the  thermal 
conversion  efficiencies  have  been  significantly  improved,  the 
average  coolant  inlet  temperature  contribution  rate  to  the  output 
power  CRcooiant  inlet  temperature  is  -0.77  W/I<  means  the  output 
power  increases  when  the  coolant  inlet  temperature  decreases,  and 
vice  versa),  whose  absolute  value  is  higher  than  CRexhaust  inlet 

temperature- 

Although  reducing  the  inlet  temperature  of  the  coolant  is  sig¬ 
nificant  to  improve  the  output  performance  of  TEG,  since  the  TEG 


Fig.  12.  Variations  of  the  output  power,  the  TEM  conversion  efficiency  and  the  TEG 
conversion  efficiency  for  the  whole  TEG  model  with  the  inlet  temperature  of  exhaust. 


Fig.  13.  Variations  of  the  output  power,  the  TEM  conversion  efficiency  and  the  TEG 
conversion  efficiency  for  the  whole  TEG  model  with  the  inlet  temperature  of  coolant. 


cooling  cycle  uses  a  branch  of  the  engine  cooling  system,  adding  a 
separate  cooling  system  or  an  intercooled  equipment  for  the 
existing  TEG  will  result  in  additional  power  consumption,  and  it 
needs  to  re-analyse  the  contribution  of  TEG  to  the  fuel  economy  of 
vehicle. 

Fig.  14  shows  the  variations  of  the  output  power  with  the 
exhaust  inlet  mass  flow  rate  under  different  exhaust  inlet  tem¬ 
peratures,  the  thickness  of  the  insulation  layer  is  10  mm  and  b  is 
15  mm.  As  the  exhaust  mass  flow  rate  increases,  the  output  power 
of  TEG  increases  obviously.  The  higher  the  exhaust  inlet  temper¬ 
ature  is,  the  faster  the  output  power  rises.  The  variations  of  the 
TEM  and  the  TEG  conversion  efficiencies  with  the  exhaust  inlet 
mass  flow  rate  under  different  exhaust  inlet  temperatures  are 
shown  in  Fig.  15.  Similar  to  Fig.  14,  the  TEM  and  the  TEG  conver¬ 
sion  efficiencies  increase  with  the  increasing  of  the  exhaust  inlet 
mass  flow  rate.  In  Fig.  11,  the  average  exhaust  inlet  mass  flow 
contribution  rate  to  the  output  power  CRexhaust  inlet  mass  flow  rate  is 
985  W/(kg/s),  1076.1  W/(kg/s)  and  1170.8  W/(kg/s)  respectively 
corresponding  to  the  exhaust  inlet  temperature  of  773  K,  793  K 
and  813  K.  It  shows  that  improving  the  exhaust  inlet  mass  flow 
rate  can  obviously  increase  the  output  power  under  higher 
exhaust  inlet  temperatures. 


Fig.  14.  Variations  of  the  TEG  output  power  with  the  exhaust  inlet  mass  flow  rate 
under  different  exhaust  inlet  temperatures. 
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Fig.  15.  Variations  of  the  TEM  and  the  TEG  conversion  efficiencies  with  the  exhaust  Fig.  17.  Variations  of  the  TEM  and  the  TEG  conversion  efficiencies  with  the  coolant 
inlet  mass  flow  rate  under  different  exhaust  inlet  temperatures.  inlet  mass  flow  rate  under  different  coolant  inlet  temperatures. 


As  is  shown  in  Fig.  8,  due  to  the  relatively  high  specific  heat 
capacity  of  the  coolant,  the  temperature  variation  of  coolant  is  less 
than  1  K.  Fig.  16  shows  the  variations  of  the  TEG  output  power  with 
the  coolant  inlet  mass  flow  rate  under  different  coolant  inlet  tem¬ 
peratures,  and  Fig.  17  indicates  the  variations  of  the  TEM  and  the 
TEG  conversion  efficiencies  with  the  coolant  inlet  mass  flow  rate 
under  different  coolant  inlet  temperatures,  the  thickness  of  the 
insulation  layer  is  10  mm  and  b  is  15  mm  for  both  of  the  two  figures. 
As  it  can  be  seen,  either  for  the  output  power,  or  the  conversion 
efficiencies,  the  influence  of  the  coolant  inlet  mass  flow  rate  is  very 
small  after  about  0.165  kg/s  (the  velocity  of  coolant  is  0.5  m/s), 
which  is  due  to  the  relatively  high  specific  heat  capacity  of  the 
coolant.  Before  this  value,  the  output  power  and  the  conversion 
efficiencies  change  obviously  with  the  increasing  of  the  coolant 
inlet  mass  flow  rate,  which  is  called  the  “sensitive  range”.  The 
average  coolant  inlet  mass  flow  contribution  rate  to  the  output 
power  CRcooiant  inlet  mass  flow  rate  in  the  sensitive  range  is  41.3  W/(kg/ 
s),  46.9  W/(kg/s)  and  52.8  W/(kg/s)  respectively  corresponding  to 
the  coolant  inlet  temperature  of  353  K,  333  K  and  313  K  (certainly, 
since  the  density  of  exhaust  and  the  density  of  coolant  is  not  in  an 


order  of  magnitude,  this  contribution  rate  should  not  be  compared 
with  CRexhaust  inlet  mass  flow  rate)- 

Based  on  the  above  analyses,  if  maintaining  the  existing 
equipment  and  structure,  increasing  the  exhaust  inlet  tempera¬ 
ture  and  the  mass  flow  rate,  as  well  as  lowering  the  coolant  inlet 
temperature  are  useful  methods  to  improve  the  output  perfor¬ 
mance  of  TEG.  The  method  of  increasing  the  exhaust  inlet  tem¬ 
perature  may  be  implemented  by  laying  TEG  on  the  front 
position  of  exhaust  tube  as  possible,  laying  an  insulation  layer 
outside  the  shell  of  TEG,  etc.  All  these  methods  will  not  generate 
any  additional  engine  output  power  consumptions.  Increasing 
the  exhaust  inlet  mass  flow  rate  and  reducing  the  coolant  inlet 
temperature  both  need  to  consume  extra  engine  output  power. 
The  net  output  power  of  TEG  Pne t  should  be  equal  to  the  output 
power  generated  by  TEG  minus  the  extra  output  power  con¬ 
sumption,  so  the  appropriate  exhaust  inlet  mass  flow  rate  and 
the  coolant  inlet  temperature  should  be  chosen  in  reason  to 
maximize  the  Pnet. 

Fig.  18  shows  the  variations  of  the  output  power,  the  TEM  con¬ 
version  efficiency  and  the  TEG  conversion  efficiency  with  b ,  which 
is  the  distance  from  the  base  of  the  heat  transfer  fin  to  TEM,  the 
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Fig.  16.  Variations  of  the  TEG  output  power  with  the  coolant  inlet  mass  flow  rate  Fig.  18.  Variations  of  the  output  power,  the  TEM  conversion  efficiency  and  the  TEG 
under  different  coolant  inlet  temperatures.  conversion  efficiency  with  b. 
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Table  2 

Summary  of  the  different  CRvar- 


CRinsulation  layer  thickness 
CRexhaust  inlet  temperature 
CRcoolant  inlet  temperature 
CRexhaust  inlet  mass  flow 

0.49  W/mm 
0.52  W/K 
-0.77  W/I< 

985  W/(kg/s) 

1076.1  W/(kg/s) 

1170.8  W/(kg/s) 

rate(773  K,  793  K,  813  I<) 

CRcoolant  inlet  mass  flow 

41.3  W/(kg/s) 

46.9  W/(kg/s) 

52.8  W/(kg/s) 

rate(353  K,  333  K,  313  K 
before  0.165  kg/s) 

CRb 

-0.24  W/mm 

thickness  of  the  insulation  layer  is  10  mm.  The  output  power  and 
the  conversion  efficiencies  decrease  to  some  extent  with  the 
increasing  of  b.  Since  the  radial  height  of  the  fin  in  this  numerical 
model  is  limited,  the  variation  range  of  b  is  only  within  a  small 
interval.  If  the  fin  is  higher,  this  effect  will  be  more  apparent.  The 
average  b  contribution  rate  to  the  output  power  CRb  in  this  simu¬ 
lation  is  -0.24  W/mm.  It  can  be  concluded  that  in  the  case  of  an 
engineering  perspective,  TEM  should  be  arranged  in  the  position 
which  is  relatively  close  to  the  roots  of  the  fin  as  possible  to  achieve 
better  output  performance. 

Table  2  is  the  summary  of  different  CRvar  obtained  from  the 
above  analyses. 

5.  Conclusion 

Based  on  the  thermoelectric  power  generation  theory,  non¬ 
equilibrium  thermodynamics  theory,  heat  transfer  theory  and 
finite  difference  principle,  a  numerical  model  for  TEG  with  cy¬ 
lindrical  shell  and  straight  fins  is  established.  According  to  the 
different  structures  of  exhaust  tube,  cylindrical  shell  and  heat 
transfer  fins,  grids  are  divided  separately  in  order  to  analyse  the 
numerical  model,  and  the  contact  effects  are  also  taken  into 
consideration.  Along  the  axial  direction,  the  temperatures  of  the 
exhaust  and  the  coolant  change  linearly.  Due  to  the  relatively 
high  specific  heat  capacity  of  the  coolant,  the  temperature  vari¬ 
ation  of  the  coolant  is  negligible.  Laying  an  insulation  layer 
outside  the  TEG  cylindrical  shell  can  effectively  improve  the 
output  power  of  TEG.  The  influences  of  the  inlet  temperature  and 
mass  flow  rate  of  the  exhaust  and  the  coolant  are  researched 
separately.  By  comparing  the  different  variable  contribution  rates 
to  the  output  power,  it  proves  that  increasing  the  inlet  temper¬ 
ature  and  the  inlet  mass  flow  rate  of  the  exhaust,  as  well  as 
reducing  the  inlet  temperature  of  the  coolant  are  useful  methods 
to  improve  the  performance  of  TEG  in  a  certain  extent.  However, 
some  methods  make  additional  power  losses,  and  the  input 
parameters  need  to  be  chosen  appropriately  to  maximize  the 
net  output  power  of  TEG.  In  addition,  if  conditions  permit, 
TEM  should  be  arranged  near  the  root  of  heat  transfer  fin  as 
possible. 

The  heat  transfer  numerical  model  established  in  this  paper  can 
not  only  be  used  for  the  newly  designed  TEG  with  cylindrical  shell 
and  straight  fins,  but  also  applied  to  the  other  types  of  TEG.  The 
versatility  of  the  numerical  model  can  evaluate  the  output  perfor¬ 
mance  of  TEG  on  the  early  stages  of  the  designing  process,  and 
optimize  the  structural  parameters,  the  external  boundary  condi¬ 
tions  and  the  related  input  parameters. 

The  next  step  is  to  establish  a  more  accurate  model  for  TEG  with 
cylindrical  shell  and  straight  fins  in  CFD  simulation  software 
environment,  and  make  some  research  on  the  heat  exchange  pro¬ 
cess.  Constructing  and  operating  a  scaled  real  TEG  with  cylindrical 
shell  and  straight  fins,  as  well  as  testing  the  performance  of  the  TEG 
to  validate  the  reliability  of  the  numerical  model  are  also  included 


in  the  program.  By  using  the  engine  simulation  software,  the 
transient  output  data  of  exhaust  can  be  obtained,  so  the  real-time 
simulations  of  TEG  will  be  carried  out  to  evaluate  its  contribution 
to  the  fuel  economy  of  vehicle. 

Nomenclature 

^TEieg  cross  sectional  area  of  thermocouple  leg  (m2) 
a  side  length  of  TEM  (m) 

b  radial  distance  from  the  root  of  the  heat  transfer  fin  to 

TEM  (m) 

Cp  constant-pressure  specific  heat  (J/(kg  I<)) 

CR  the  average  variable  contribution  rate  to  the  output 
power 

d  diameter  (m) 

Gra ir  Grashof  number  of  air 

h  heat  convection  coefficient  (W/(m2  K)) 

/  current  (A) 

I<  thermal  conduction  (W/m) 

L  total  axial  length  of  fin  (m) 

Lrmeg  length  of  thermocouple  leg  (m) 
m  mass  flow  rate  (kg/s) 

Nfin  number  of  fins 

Nuf  Nusselt  number  of  fluid 

P  output  power  (W) 

Prf  Prandtl  number  of  fluid 

Qh  absorbing  heat  of  TEM  ( W) 

Qc  releasing  heat  of  TEM  (W) 

q  heat  mass  flow  (W) 

R  internal  resistance  of  thermocouple  (Q) 

Ri  load  resistance  (Q) 

^copper  resistance  of  copper  stripe  (Q) 

Ref  Reynolds  number  of  fluid 

T  absolute  temperature  (K) 

t  thickness  (m) 

a  Seebeck  coefficient  (V/I<) 

X  thermal  conductivity  (W/(m  I<)) 

p  density  (kg/m3) 

j'air  kinematic  viscosity  of  air  (m2/s) 

7]  heat  conversion  efficiency  (W) 

a  electrical  resistivity  (Q  m) 

Ax  axial  length  of  the  control  volume  (m) 

Ay  radial  length  of  the  control  volume  on  fins  (m) 

Subscripts 

cf  cooling  fluid 

cf-wall  between  the  cooling  fluid  and  the  tube  wall 

con_ex  convection  of  the  exhaust 

con_cf  convection  of  the  cooling  fluid 

ct  cooling  tube 

cw  cooling  tube  wall 

cyl  cylinder  shell 

ex  exhaust 

ex_tube  exhaust  tube 

ex-wall  between  the  exhaust  and  the  tube  wall 
fin_air  between  the  fin  and  air 

i  number  of  axial  control  volume 

ins  insulation  layer 

k  number  of  radial  control  volume  for  ith  CV 

surf  surface 

tc  cold  side  of  TEM 

th  hot  side  of  TEM 

tw  tube  wall 

var  variable 
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